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We investigate superfluidity in the interior-gap states proposed by Liu and Wilczek. At weak 
coupling, we find the gapless interior-gap state unstable in physically accessible regimes of the 
parameter space, where the superfluid density is shown to be always negative. We therefore conclude 
that the spatially-uniform interior-gap phase is extremely unstable unless it is fully gapped; in this 
case, however, the state is rather similar to conventional BCS states. 
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I. INTRODUCTION 

The advances of techniques in manipulating dilute 
cold atoms in traps have opened up many fascinat- 
ing new possibilities to condensed matter physics in 
recent years, p] The Bose-Einstein condensation of di- 
lute alkali atoms has provided new Bose-condensed sys- 
tems previously known only for liquid 4 He. Experi- 
mentalists are now investigating the possibility of re- 
alizing the Bardeen-Cooper-Schrieffer (BCS) states in 
trapped Fermionic atoms. 0, 0, Q In the simplest case of 
spin-singlet BCS states, as in ordinary superconductors 
such as aluminum, the "conventional" Cooper pairs are 
formed from Fermions with opposite spins. In atom 
traps, a corresponding realization could be two species 
of Fermions which pair via the inter-species interaction. 
The two species of Fermions can be the same kind of 
atoms (say, 6 Li) in different hyperfine states, or, in the 
most general case, different Fermionic atoms. 

Conventionally, the BCS state is a condensate of 
Cooper pairs consisting of Fermions with equal mass 
and opposite spins from states with a single Fermi sur- 
face. In considering superconductivity in systems con- 
taining fcrromagnctically coupled paramagnetic impuri- 
ties, Fulde and Ferrell, and independently Larkin and 
Ovchinnikov0 (FFLO), studied Cooper-pairing in sys- 
tems with mismatched Fermi surfaces. They found a spa- 
tially non-uniform phase can be more favored than the 
BCS states in certain ranges of temperature and Fermi- 
surface mismatches. Recently Liu and Wilczek (LW), 
motivated by recent developments in atomic physics, 
studied a system of interacting Fermions consisting of 
two species of particles with unequal masses and mis- 
matched Fermi surfaces. Assuming pairing only in the 
vicinity of the smaller Fermi surface, they claimed that at 
weak coupling there exists a novel superconducting state, 
which is spatially uniform, could be energetically more 
favorable than the FFLO state if the coupling strength 
is above a critical value. This state is characterized by 
the coexistence of superfluid and normal components at 
zero temperature that are separated by a momentum gap 
between the two Fermi surfaces (thus the "interior-gap" 
state). LW concluded that the critical coupling strength 
for the interior-gap state could be vanishingly small when 



the mass difference between the two Fermion species is 
large. 

In this paper, we re-examine the possibility of an 
interior-gap phase having in mind cold dilute Fermionic 
atoms in traps. Instead of postulating pairing only 
around the smaller Fermi surfaces, [(| here we consider 
a short-range attractive interaction arising from the low 
energy scattering between Fermions of unlike species. We 
derive and solve the gap equation, and then further ex- 
amine the current response of the interior-gap state in 
the weak coupling limit. Surprisingly, we find the state 
unstable in all physically accessible regimes unless it is 
fully gapped. 



II. THE MODEL AND THE GAP EQUATION 

We consider an atom trap which contains two species 
of Fermions with mass rrih > mi. We assume that the 
gas is dilute, with the Fermions interacting via a short- 
range interaction (i.e., both the scattering length and the 
range of the interaction between the two Fermion species 
are much less than the interparticle distances). This in- 
teraction can therefore be modeled by a delta-function 
like coupling of strength —g > which ultimately causes 
pairing between the heavy (/i) and the light (I) particles. 
The system is thus described by the Hamiltonian 

H = Yl £p" a P" a P« + f J2 a l>h a -pn a -pi a ph- (1) 

pa=h : l P,p' 

Here a pa and a pa annihilates and creates an a species 
particle with momentum p. For simplicity, we ignore the 
trapping potential and take the parabolic dispersions for 
the particles £ pa = (p 2 — p 2 Fa )/2m a with pFa the Fermi 
momentum. 

Since the gas is dilute, we shall make mean-field ap- 
proximation and define the space-independent order pa- 
rameter A = — gJ2 P ( a ph a -pi) which is chosen to be real. 
The quasiparticle dispersions are then obtained from the 
standard procedure, |5j yielding 



E phJ = ±^L_^L + ^ ^U±ipl j 2 + A2 . (2) 
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FIG. 1: The right hand side of the gap equation [in units of 
2ir I (pFhm r )] for mj, = mj; (a) shows the three-dimensional 
plot, and (b) the profiles at p F i/p F h=0-5 (dot-dashed line), 1.0 
(full line), and 1.5 (dashed line). The ridges in (a) correspond 
to the critical value A c in Eq. (|1J which appear as peaks in 
(b). 



The gap equation is derived from the Hamiltonian upon 
minimizing the free energy with respect to A. However, 
due to the short-range nature of the pairing interaction, 
the resulting formula has ultraviolet divergence. To re- 
solve the problem, we eliminate the coupling constant g 
in favor of the scattering length a for two-particle scat- 
tering (between h and I) in vacuum. The regularized gap 
equation readsQ 



m r ^ — ^ 



1 



E h + E, 



(1 - f(E h ) - f(Ei)) - -q 



1 



,(3) 



where m r = m,hmij(mh + mi) is the reduced mass, a 
the s-wave scattering length; / is the Fermi function and 
= p 2 /2m a the particle dispersions in vacuum. Here 
and in below, for brevity, we omit the subscript p when 
no confusion would likely arise. 

We solve the gap equation at zero temperature by plot- 
ting the right hand side (RHS) of Eq. © for different val- 
ues of mi I rah over the parameter space (pp,A), where 
pp = pFi/PFh and A = |A|/2|/j,^| with [ih the chemical 
potential of the heavy particles. Since a < for attrac- 
tive interactions (without formation of h-l bound states 




FIG. 2: The quasiparticle dispersions Eh and — Ei for mi = 
7mi and pFh = 2-Opfi', the normal-state dispersions £h and 
— are also shown for the heavy (dotted line) and the light 
(dashed line) species. Here A = 0.6|/ih| < A c = 0.99|^ih|. It 
is seen that Eh intersects the p-axis at the points pi and p2, 
which correspond to the regions of gapless excitations. 



in vacuum), we seek solutions of A such that the RHS is 
positive. FigureHshows typical plots for the RHS of the 
gap equation with rrih = mi. The profile at pfi/pfh—^- m 
Fig.^6) thus corresponds to the usual BCS case. In this 
case the RHS diverges as A — > 0, corresponding to the 
well-known fact that an infmitesimally small attractive 
interaction can lead to pairing. However, two differences 
from the usual BCS states occur here. Firstly, the RHS 
becomes negative at large A, indicating no solution for 
attractive interaction. Secondly, when PFh ^ Pfu as dis- 
played in Fig.^6), a finite peak (maximum as a function 
of A) arises in the RHS of the gap equation at finite A. 
In the three-dimensional plots of Fig. |IJa), these peaks 
appear as ridges over the (pf ; A) plane. The existence 
of these peaks (ridges) implies that the system requires a 
critical coupling strength for solutions to the gap equa- 
tion when there are Fermi surface mismatches. Another 
significance of these peaks and ridges is related to the 
quasiparticle properties of the system, as we shall now 
explain. 

An important feature of the interior-gap state as pro- 
posed by LW is the coexistence of the superfluid and the 
normal components at zero temperature; in other words, 
there exists gapless excitations. However, as one can 
check from Eq. J5J), this is possible only when the mag- 
nitude of the order parameter |A| is smaller than 



\PFh-P 2 Fl\ 



(4) 



A^/m h mi 

When |A| < A c , depending on the relative magnitude of 
PFh and pfu either Eh or Ei crosses zero at the points 
(see Fig. EJ 

1/2 



Pi, 2 



PFh + Pfi T ViPFh - Pfi) 2 - Wm h miA 2 



(5) 

In the event |A| > A c , Eh,i both stay positive for all val- 
ues of p. There are then no gapless excitations and the 
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interior-gap state has only the superfluid component at 
zero temperature. The peaks and ridges in the plots for 
the RHS turn out to locate exactly at the value |A| = A c 
of Eq. Q. Indeed since |A| = A c draws the boundary 
between the gapped and the gapless regions over the pa- 
rameter space, one can expect a qualitative change when 
plotting the RHS across this line. We will be interested 
mainly in the gapless regions, namely where |A| < A c . 



III. CURRENT RESPONSE 

When |A| < A c the interior-gap state has gapless ex- 
citations from the quasiparticle branch whose dispersion 
crosses zero. In the FFLO ground state, the quasipar- 
ticles produce a flow which cancels exactly the current 
due to the finite-momentum Cooper pairs. This lowers 
the free energy and stabilizes the FFLO states. Q In the 
interior-gap state, the Cooper pairs are stationary in the 
absence of any superflow. It is therefore of interest to in- 
vestigate the current response of the interior gap state; in 
particular, the effects of the quasiparticlcs in the current- 
carrying states. 

In the static case the quasiparticle distributions of the 
interior-gap phase are given by the expression 

n° pa = u 2 p f(E pa ) + v 2 p f(-E p& ) , (6) 

where a is the species other than a and the coherence 
factors are 



2 E h Ei+^h 2 
U„ = = , V- = 1 

p E h + Ei E h + Ei' p 



(7) 



In the presence of a small superfluid velocity w, the quasi- 
particle energies are shifted by +(p • w) so that the quasi- 
particle distribution functions of Eq. © become 



i p « = u 2 p f{E pa + p • w) + vlf(-E pS + p • 



(8) 



The number current can be decomposed as usual into two 
parts [f| 



(9) 
(10) 



where the superscripts p and d indicate the paramag- 
netic and the diamagnetic components; p^ d are the cor- 
responding number densities. In the following we shall 
examine Eqs. © and IjlUI) in the limit of small w at zero 
temperature. 

Let us consider the case when pph > Pfi and | A| < A c . 
Then, as can be checked from Eh < for p\ < 
p < P2 while Ei > always. Therefore the I branch 
of quasiparticles is always empty, i.e., f{E{] — and 
f(—Ei) = 1, at zero temperature. In the presence of a 



small superfluid velocity w, since Ei never changes sign, 
one has 



f(±E l +p-w)-f(±E l ) = 0. 



(11) 



Using the fact that 3 p a = when w = 0, we obtain for 
small superfluid velocity the paramagnetic currents 

J « = ^Ep(v-«p a ) 

~ — V p 4>l [f(r) a E h + p ■ w) - f( Va E h )}(12) 
p 

where for the heavy and the light branches 



and rj a — 



for a — 



I. 



(13) 



In arriving at the final expression in Eq. I|12f) . we have 
used 111 111 for small w. 

To leading order in the superfluid velocity, the term 
in the square brackets in Eq. i|12|) is proportional to the 
delta function 5(Eh), which vanishes except at p = p± 
and p — pi- Replacing the sum by integral, one thus 
obtains 



ptD 



2<P a 2) 



(14) 



where Di = l/\dEh/dp\i and the subscripts i = 1,2 in- 
dicate evaluating at the points p±, p2- Note that p v a is 
always negative, which means that the paramagnetic cur- 
rent always flows in the opposite direction to w. 

To leading order in w, the diamagnetic current is sim- 
ply 



p 



pa 



w . 



(15) 



which is always in the same direction as w. When pph > 
Pfi the diamagnetic number-densities can be expressed 
as 



Ph - 1 - 6^ + 2^ 



2„,2 



dpp u 



dppV . (16) 



Note that p d a is always positive since the diamagnetic 
current J^, is always in the same direction as w. 

Similar analysis as above can also be done for the case 
of pFh < Pfi- As it turns out, we find that the total su- 
perfluid number densities p™ = p^+Pa f° r l ne two species 
are always identical (see Appendix). In other words, 
despite the seemingly asymmetric properties of the two 
species, the superflow turns out to be quite conventional 
- the particles flow together in pairs. The superfluid mass 
density is thus simply p M = m h Ph+ m iPi = {mh+mi)Ph- 

As noted earlier, p d a is always positive while fP a is 
always negative; as a result, the sign of pu is deter- 
mined from the competition between the paramagnetic 
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FIG. 3: Summary of our results for the interior-gap states 
for mi/rrih = (a) 1, (b) 1/7, and (c) 1/100 over the (pF,A) 
plane. The grey areas are regions where the gap equation has 
solutions for attractive interactions. The dark areas mark the 
regions where the superfluid mass density is positive in the 
gapless regions. The dash-dotted line depicts the critical line 
determined from Eq. (|IJ, which also corresponds to the ridges 
in the RHS of the gap equation (see Fig. 0[a)). 



and the diamagnetic components. Quite unexpectedly, 
we find the superfluid mass density of the interior-gap 
phase stays negative in almost all regions of the param- 
eter space. This can also be confirmed analytically for 
small A. A negative superfluid mass density implies 
that the interior-gap phase is unstable towards a spon- 
taneously generated superfluid velocity or a phase gra- 
dient, since the kinetic energy is \pm~w 2 < 0. Our re- 
sults for mi/nih —I, 1/7, and 1/100 are shown in Fig.|3] 
From Fig.[3a) one observes that for to^ = to/ the regions 
where pu > and the regions where the gap equation 
has solutions are completely disjoint to each other in the 



FIG. 4: The RHS of the gap equation [in units of 2ir / (pFim r )] 
for nrih = 7m; at Pfi = 2pfh- The dashed line indicates the 
value of A below which the superfluid mass density is positive 
(cf. Fig. To the right of the peak, the state becomes 

fully gapped, so that the mass density is again positive triv- 
ially. 



gapless region. Namely for = to/, except in the fully 
gapped region, there is no where in the parameter space 
(pF, A) the interior-gap state being stable. This is con- 
sistent with the FFLO results. For larger mass ratio 
at to/j = 7to/, as shown in Fig. &b), there are regions 
in the gapless area where pM stays positive and the gap 
equation has solutions. To estimate the coupling strength 
required, we plot in Fig. 21 the profile of the RHS of the 
gap equation at ppi = %PFh- Disregarding the fact that 
there are two solutions to the gap equation in this case 
(see next Section for a discussion), we find the coupling- 
strength necessary for pm > around |pir/a| ~ 1.9. At 
such strong couplings, however, the validity of the present 
mean-field calculation is doubtful. For even larger mass 
ratio at rrih — 100to/, the situation stays much the same 
(see Fig. EUc)); at ppi = ^Pfk we find \pf\o\ ~ 0.7 for 
Pm to be positive. Even though the coupling strength 
does go down with increasing mass ratio, it is unrealistic 
to consider even higher values of to^/to/. Note also that 
the regions of pu > occurs for ppi > PFh, instead of 
Pfi < PFh as suggested by LW.[|| 

From the sign of the superfluid density, our re- 
sults show that in the interior-gap phase the paramag- 
netic component overwhelms the diamagnetic component 
whenever there are gapless excitations. It is only when 
the quasiparticles are gapped out the number density in 
the interior-gap state would stay positive at weak cou- 
plings. 



IV. SUMMARY AND DISCUSSION 

Within the range of validity of our calculation, we find 
the interior-gap phase unstable in physically accessible 
regions of the parameter space; the superfluid mass den- 
sity is always negative at weak couplings, rendering a 
uniform phase unstable. Therefore, starting from the 
normal state, upon increasing coupling strength the sys- 
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tern considered here would/may first enter an inhomoge- 
neous FFLO state which persists until the onset of a fully 
gapped "interior-gap" state; a uniform gapless interior- 
gap state has no room in the intervening regimes. Indeed 
this may have been expected from continuity: As was no- 
ticed, due to the existence of the peaks and ridges in the 
RHS, there can be two solutions to the gap equation (see 
Fig.0. Since at pfk = Pfi the only solution is always in 
the gapped region, from continuity one expects the sta- 
ble solutions when pph ^ Pfi to be those above the V 
shape boundaries in Fig. In other words, continuity 
implies that it is always the gapful solutions that are sta- 
ble. This includes the large \pfi<A regimes discussed near 
the end of the last section. For this fully gapped phase, 
there are no quasiparticles at zero temperature and the 
number densities of the h and the I branches are always 
equal. [Tol] This state thus behaves just like an ordinary 
BCS state. 

This research was supported by NSC of Taiwan under 
grant number NSC 91-2112-M-001-063. 
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Noting that Eh(pi,2) = 0, one can derive from @ 
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Therefore (|A 2|) becomes 
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(A 2) 
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For the diamagnetic number densities, at small super- 
fluid velocity the difference between the two species is 
simply that of their quasiparticle occupation numbers. 
As one can find trivially from Eq. 



Appendix 

In this appendix we show that the superfluid number 
densities for the heavy and the light species are identical 
at small superfluid velocities. 

Let us consider first the case of pfk > Pfi- We note 
first that from Eqs. 101 and 
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Combining Eqs. (|A 4(1 and 1\A 5|l . we conclude that p\ — 
p'l = 0, namely that the superfluid number densities of 
the two species are identical. 

For pph < Pfi similar calculations as above lead to 
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Substituting the above expression into Eq. I|14[l . we ob- Therefore, again, one finds p]l — in this case. 
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